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Abstract 

The dynamics and radiation of positronium is investigated in intense laser fields. Our two-body 
quantum mechanical treatment displays the tunneling, free-evolution and recollision dynamics of 
electron and positron both in the oscillating laser electric and laser magnetic field components. 
In spite of significant momentum transfer of the numerous incoming laser photons, recollisions of 
both particles are shown to occur automatically after tunneling ionization, along with substantial 
x-ray and gamma-ray emission during recombination and annihilation processes. 
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The highly nonlinear interaction of gaseous atoms with super intense laser pulses has been 
demonstrated to give rise to the emission of coherent high-harmonic generation (HHG) 
up to the x-ray regime, with numerous applications in high-resolution spectroscopy and 
diagnostics (see e.g. Since the kinetic energy and consequently the frequency of the 

emitted radiation of a particle increases with rising laser intensity, considerable interest has 
been directed towards understanding the complex relativistic quantum dynamics of atoms 
and ions in ultra intense laser pulses 2J. Once electrons reach velocities nonnegligible to that 
of light, however, the magnetically induced Lorentz force or in other words the momentum 
transfer of the numerous incoming photons induce a separation of electrons and ionic core 
in the laser propagation direction, resulting in a strong reduction of high harmonic yields 
While highly charged ions and crystals were studied towards a reduction of this 
effect p , there is still a clear lack for an efficient system, where radiation pressure does not 
induce substantial ionization in the laser propagation direction and thus reduce coherent 
high frequency generation. On a different front there is also a quest for physics beyond 
atomic and c assical plasma physics in ultra intense lasers, such as nuclear reactions Jf| and 
QED effects J. 

In this letter the dynamics and high harmonic generation of positronium is investigated in 
intense laser fields. The two-body system is shown to display unique properties: While 
tunnel-ionization of electron and positron may occur almost oppositely in the laser polar- 
isation direction, both particles sense the identical drift in the laser propagation direction 
due to their equal magnitudes of mass and charge (see Fig. [TJ. Periodic electron-positron 
recollisions are shown to occur automatically head-on in spite of the influence of the Lorentz 
force. In addition to substantial coherent x-ray generation during periodic electron-positron 
recombinations we predict gamma radiation in the less likely events of laser-enhanced anni- 
hilations of both particles. 

Positronium consists of an electron and a positron and is known to be unstable. While ortho- 
positronium annihilates into three photons with a lifetime of 1.4 ■ 10~ 7 s, para-positronium 
does so with two photons and a lifetime of 1.25 • 10^ 10 s jsj]. The presence of a strong 
laser field may induce substantial reductions l| or enhancements f]~oL fll| of the annihilation 
process into gamma photons. However, even for the shorter lifetime of para-positronium it is 
sufficient for the interaction with many cycles of a femto second laser pulse with interaction 

n 

lengths not exceeding much the centimeter range [10]. 
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The laser field shall be described by the vector potential A(x, t) propagating in the z direc- 
tion and linearly polarized along the x direction. Being interested in the tunneling regime 
with moderately intense laser field strengths we restrict ourselves to Schrodinger dynamics 
in two dimensions, though beyond the dipole approximation. The quantum dynamics of 
positronium in the laser field, taking fully into account of the Lorentz force due to the laser 
magnetic field component, is thus governed by the following Hamiltonian (in atomic units 
as throughout the article): 



(i) 
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where, x$ and — iV Xj , (i £ { e ,p})> are the operators of the coordinates and momenta of the 
electron and positron, respectively. Further we introduce relative r = (x, y, z) = x e — x p and 
center of mass R = (X, Y, Z) = Xe + Xp coordinates. In accordance with the Ehrenfest theo- 
rem, the center of mass transversal canonical momenta —idx, —idy as well as i (d t + cdz) 
turn out to be conserved quantities, i.e. commute with H. The eigenvalues are defined to be 
P x ,P y ,g, respectively, and £ is a separation constant that can be understood as the energy 
of the system before the interaction. Therefore, we consider the ansatz 

$(r, R, r) = exp (i (p x X + P y Y - 9 -Z - 5r)) 0(r, r) (2) 

singling out eigenfunctions of the conserved quantities in the wave function and introducing 
the running time r = t — —. Employing a 1/c expansion of the vector potential as function 
of the coordinates of the relative motion of the particles, we obtain the following equation 
for d>: 



i(l-B z )d T <f)(r,T) 
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Here, the center of mass velocity is introduced via V = P/2 with B = V/c and 
P = (P x ,P y ,(£ — g)/c). The applied expansion takes into account of the leading multi- 
pole operators describing the magnetically induced relative motion of electron and positron 
in the laser propagation direction. Since P x is a conserved quantity, this operator describes 



an oscillation of the relative coordinate in the propagation direction while for atoms the 
analogous term would lead to a drift. For consistency, the higher-order term for the center- 
of-mass motion —■Asd^. in Eq. (jSJ) will be neglected in the following. 

We proceed by carrying out a transformation to the length gauge: exp(— ix(r, r))0 = \I/ 
with x(r, r) = A(r) ■ r/c. In this gauge we note the full velocity dependence of the transition 
matrix elements to be discussed later. The resulting time dependent Schrodinger equation 
reads 



i(l-B z )d r *(r,T) = < -V 2 r + I p 



-E(t)(1-B z )x-B x E(t)z-±\*( 



(4) 



with electric field E(r) = — \d T A and positronium ionization potential I p . Further we 
employ the reasonable assumptions proposed earlier by Lewenstein and co-workers Q| for 
atomic HHG in the dipole approximation, being: a) The contribution of all bound states 
except the ground state is neglected; b) The depletion of the ground state is neglected; c) 
In the continuum the electron is treated as a free particle in the laser field. Then the time 
dependent wave function can be written as 

/oo 
d 3 p6(p,r)|p> (5) 
-oo 

Here |0) is the ground state and b(p, r) denotes the amplitude of the corresponding contin- 
uum state |p). Integration of Eq. (J3J) with the ansatz © yields for the ionization amplitude 
b: 

6(p,r) = — l — [ T ( P -A(tO/c-A(t')|^(t')|0> 

x exp ( — iS (p, t, t')) c£t', (6) 

where the quasiclassical action is S(p, r, t') = J T , dr" {(p — A(r")/c — A(r")) 2 + J p } /(l — 
B z ). Here, the new variable p = p + A(r)/c + A(r) is introduced with A(r) = A(r)e 2 = 
B x A(t)/ [c(1 — B z )] e z and Hi{t) = —E{r)x{\ — B z ) — B x E(r)z. In solving the equation 
for b we have neglected all terms involving continuum-continuum transitions of no interest 

n 

for the high-frequency part of harmonic generation as in 12 • 

We calculate the dipole moment x^r) = (\&(r)| r |\I/(r)) of which the Fourier transform 
x. d ([uj — k - V]/[l — B z ]) yields the Doppler shifted harmonic spectrum, where lj and k denote 



the emitted radiation frequency and wave vector, respectively. In the tunneling regime the 
integrals can be obtained using the saddle point method, V P S (p, r, r') = 0. Sustaining all 
momenta up to first order in 1/c yields: p x (r, r') = f T , dr" A(r") / c(t — r f ),p y = 0,p z (r, r' 
B x p x (t,t'). The matrix elements can be evaluated analytically employing the relation 
d(p) = (d x ,d y ,d z ) = (p| r |0) = Ai (I p ) 5 ^ 4 p/(vr (p 2 + I p ) 3 ) and performing the integration 
with respect to p gives 

xd* x (p(r,r')-A(r)/c-A(r)) 
xE(r')d x ( P (t,t')-A(t')/c-A(t')) 

x exp (-iS(p(r, r'), r, r')) + c.c., (7) 

with classical action reading S(p(t,t'),t,t') = f^, dr" [ (p(r, r') — A(r")/c — A(r")) 2 + 
Ip] / (1 — 5 2 ). Taking into account the pole in the r' integration finally results in 

T-r ) 2E{t ) 
xd x {p x (r,T ) - ^-^) 

3 1% 

x exp { - iS (p(r, r ), r, r ) } + c.c. (8) 

The birth times tq are determined, for given r, by the condition p x (r, tq) = A(tq)/c. The 
leading terms in 1/c have been sustained only in the phases because they play a far less 
substantial role in the prefactor being identical to the dipole case. 

The dipole moment in Eq. (jHJ) displays essential differences in comparison with atomic 
HHG which shows especially in the real exponential factor in Eq. (jHJ): C := exp[-2/ p 3/2 (l + 
B z ) / (3£'(ro))]. We stress that in the atomic HHG case Q there is an additional contribution 
exchanging I p by + c 2 £ 4 /2 due to the drift in the laser propagation direction which can 
strongly reduce the amplitude and thus the efficiency of HHG. Here £ = \/2A/c 2 and I£ 
denotes the atomic ionization potential. The expression C, without the factor 1 + B z , 
coincides with the corresponding result within the dipole approximation for atoms. This 
means that for positronium, the magnetically-induced drift in the laser propagation direction 
has no diminishing impact on HHG. 
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The figures 121 display harmonic spectra in the laser propagation direction following the 
numerical integration of the dipole via Eq. (JJJ). In agreement with the numerical results, 
the maximal harmonic frequency 0J max is shown to critically depend on both the center-of- 
mass velocity and the observation direction via the Doppler shift with u max (l — k • V/|k|) = 
"^max^L^X ~ B z ) = 3.17U p + I p for maximal harmonic number n max and U p = c 2 £ 2 /2. 
In order to quantify the advantage of positronium for HHG in comparison to conventional 
HHG in atomic gas jets, we consider an example. In Lithium vapor being irradiated by a 
laser field at A = 0.75/xm wavelength and intensity 1 = 7- 10 16 W/cm 2 (£ = 0.17), HHG up 
to the 10 keV photon energy range can be achieved in the tunneling regime for Li 2+ with 
Ip = 122.88 eV. Magnetic field effects in this case, however, are problematic as the extent 
of the electron drift during the laser period exceeds the atomic size: A£ 2 /167r ^> ag, where 
is the Bohr radius. This substantially reduces the HHG efficiency, because of the earlier 
discussed c 2 £ 4 term being absent in positronium. Therefore, at £ = 0.12 (e.g. A ~ 50/im and 
1 = 8- 10 12 W/cm 2 ), the HHG cut-off frequency for positronium is equal to that of Li 2+ at 
the above mentioned parameters, however with a yield 15 orders of magnitude larger. Thus, 
the low densities of positronium may be compensated for by the absence of drift-reduced 
recollisions. 

Finally we discuss the possibility of 7-ray emission with positronium via laser induced an- 
nihilation. In a super strong laser field with £ > (mc 2 /fku) 1 ^ 3 « 10 2 , annihilation can occur 
via emission of n > 10 6 photons [u| along with one 7-quantum (Fig. 3a). In more moderate 
laser intensities in the tunneling regime with £ < 1 of interest here, ortho-positronium, may 
annihilate into two 7-quanta (with induced emission of an odd number of laser photons, Fig. 
3b) or spontaneously into three photons (Fig. 3c)). 

We proceed with an order-of-magnitude estimation of 7-ray emission via strong laser assisted 
annihilation of ortho-positronium as depicted in Fig. 3b in comparison to spontaneous 3- 
photon annihilation (Fig. 3c). The 2- vortex diagram in Fig. 3b, at first with bare particles, 
can be evaluated as product of the Thomsonian cross-section a ~ i\r\ with r$ being the 
classical electron radius [si] and the particle flux to yield the annihilation probability W ~ 
a ■ pc, with electron and positron density p = a^ 3 . The annihilation process is considered in 
a strong laser field here, with laser frequency being low in comparison to the characteristic 
frequency uj c = c 2 . In this situation the probability for bare particles may be multiplied by a 
factor F s which takes into account of s laser photons being emitted during the annihilation 
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process (analogous to Low theorem The laser- induced differential probability of 7-ray 

emission at frequency uj c being accompanied by the emission of s laser photons can then be 
written as: dW in d ~ 7rrlpc(du>/ Au> D )(dfl/4:7i) F s . Here we assumed that the 7-emission is 
distributed within a Atc solid angle Q and that the Doppler width Auj£, describes the lea ding 
spectral broadening. Further F s = Al(sa/3) + £ 2 (2w — l)[A 2 (saj3) — A Q (saf3)A 2 {saf3)} [11| 
is represented via generalized Bessel functions A (saf3), Ai(saf3), A 2 {sa[3) 1 a = k'A/(kq)c, 
(3 = -kk' A 2 /(8c 2 (kq) 2 ), u = (kk') 2 /(A(kq)(kq)), where = p^ - k^A 2 / \Akp)c 2 , q^ = 
p'^ — k^A 2 1 \Akp)c 2 for {/i} G {0,1,2,3} and k,k',p,p denote the laser photon, 7-quanta, 
electron and positron 4-momenta, respectively A is the 4-vector-potential. In a weak laser 
field with ^ C 1 our formula for dWi n( i matches the known result in Q|. With a similar 
heuristic argumentation we may estimate the spontaneous 3-photon annihilation probability 
without laser field at u c depicted in Fig. 3c, to give dW sp ~ aTrr^pc^du /u c )(dVL/AT[) with 
fine structure constant a = 1/c. This result for Fig. 3c , however, is well-known from the 
literature [sj]. 

For moderately intense laser fields with £ -C 1, the main contribution to dWi n d arises from 
F\ m £ 2 , and at £ ~ 1, for particular s there are arguments with F s w 1. In this latter case, 
the induced emission dW in d can exceed the spontaneous background dW sp up to a factor of 
(w c /Awo)/a Rj 10 6 for a Doppler linewidth Aujd/uj c ~ \J\ksT /uj c ) pa 10~ 4 with Boltzmann 
constant k% and temperature T. For para-positronium, however, a similar analysis shows 
that the intensity of 7-radiation via laser-induced annihilation can not exceed that via spon- 
taneous annihilation. Regarding the 7-ray yield via laser-induced annihilation, a brilliance 
at photon energy 0.5 MeV of B ~ 10 11 photons /{s ■ mm 2 mrad 2 (0.l%bandwidth)) can be 
estimated for ortho-positronium density of 10 10 cm -3 and an interaction length of 1 cm. This 
appears clearly detectable though moderate as compared to large-scale 7-ray sources in the 
same spectral region [3]. 

Concluding, positronium was shown to provide various promising dynamical features in 
strong laser fields such as energetic head-on recollisions of its constituents, with substantial 
coherent x-ray generation and 7-ray emission. 

This work has been supported by the German science foundation (Nachwuchsgruppe within 
SFB276) and an Alexander-von-Humboldt fellowship for KZH. We benefited from help on 
numerical issues by D.Diehl and A.Staudt, and especially from a seminal discussion with D. 
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Figures 
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FIG. 1: Schematic diagram displaying positronium dynamics in an intense laser field. The bound 
system depicted by the density of its wave function may ionize in the laser field. Once free, both 
electron e and positron p could be described as classical particles. Their trajectories are shown 
by the solid lines. The electric field E accelerates both particles in opposite directions while the 
Lorentz force due to the magnetic field H leads to an identical drift in the propagation direction 
kx- Without an initial center-of-mass motion these trajectories are symmetric and thus both 
particles overlap periodically giving rise to possible recombinations and annihilations and thus 
high-frequency light emission. 
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FIG. 2: High harmonic spectra of positronium in the laser propagation direction in a laser pulse 
with vector potential A = 194.0a.u. and angular frequency oj = 0.0091a.u. (corresponding to a 
laser intensity of 8.5 x 10 12 W/cm 2 and the second harmonic of a CO2 laser with wavelenth 10 fim) 
and pulse shape parameters c\ = 1.0 and C2 = 5.0. The symbols c\ and C2 denote the number of 
cycles of sin 2 turn-on/turn-off phases and of the intermediate cycles at full amplitude, respectively. 
The initial positronium velocity is for a) V x = V z = O.OOOla.u. and for b) V x = V z = 13.7a.u. 
corresponding to 10% of the speed of light. 
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FIG. 3: Feynman diagrams displaying 7-photon emission via electron-positron annihilation by a) 
multiphoton annihilation in super strong laser fields at £ > 100 with emission of n > 10 6 laser 
photons along with one 7-quantum; b) annihilation in a laser field at £ <C 100 with emission of 
laser photons along with two 7-quanta; c) annihilation of ortho-positronium with emission of three 
photons without laser field. Bold lines correspond to the electron (positron) Volkov states, i.e. in 
the presence of the laser field, dashed lines to the emitted photons. 
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